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ABSTRACT 



This study presents the results of the aerodynamic 
design of a transonic axial-flow fan for an aircraft air- 
conditioning unit. The design is based on fan design-point 
specifications and the experimental test results for a 
model of a proposed fan— inlet. 

The fan— inlet model was tested under flow conditions 
that gave Reynolds number similarity to the design inlet 
conditions. Three-dimensional pressure probes were used to 
measure the model discharge flow conditions at eight sta- 
tions. The flow was found to lack axial symmetry. There- 
fore/ several internal modifications were made. Complete 
axial symmetry was not achieved for a reasonable loss in 
total pressure. It was necessary to circumferentially aver- 
age the velocity distribution and mass— average the total- 
pressure loss in order to provide a basis for the fan 
design. 

Concurrently another design based on uniform flow and 
model losses was made. The results were nearly identical. 



Thesis by: Ronald William Pyle and Milton Ward York entitled 

Transonic Fan For Aircraft Air-Conditioning Unit. 
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1. Introduction 



An axial— flow fan may be distinguished from an axial- 
flow compressor in that a fan is generally classified as a 
high capacity, low head, single-stage unit which operates 
at a high, constant speed. Although used for different pur- 
poses than a compressor, an axial-flow fan may be designed 
or analyzed with the available theoretical and experimental 
methods for axial-flow compressors. 

One of the principal assumptions made in the aerodynamic 
design of axial-flow compressors is that of axially symmet- 
ric flow. This condition is primarily assumed to reduce the 
general flow equations from a complex three-dimensional to a 
simplified two-dimensional system. At a given radial loca- 
tion the physical requirement of this assumption is for 
circumferentially constant flow conditions ahead of any 
blade row. 

T3ie single-stage transonic fan, with its relatively 
higher flow capacity and pressure ratio over conventional 
subsonic fans, has been successfully used in modern air- 
craft systems. However, for compact systems the transonic 
turbomachine becomes more demanding of inlet flow uniform- 
ity. Thus, proper matching of an inlet— scroll to a fan 
stage becomes more critical for a given set of dimensional 
restrictions . 
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The purpose of this paper was the aerodynamic design 
of a single-stage transonic fan for an aircraft air- 
conditioning unit. The design was to be based on a list of 
design-point specifications and also on a proposed fan in- 
let that complied with the geometrical design envelope. 

A double-scale model of the proposed fan inlet was tes- 
ted. Internal modifications were made in an attempt to 
attain uniform, axially symmetric flow at a location just 
ahead of the required rotor position. Final model test 
results were used to predict the losses through the fan in- 
let and the flow conditions ahead of the fan rotor. The 
aerodynamic fan design was based on these conditions. 

A similar fan design was made -which was based on assumed 
% 

uniform, axially symmetric flow ahead of the rotor. Th°us, 
a comparison could be made between the two designs, and 
also with the design specifications. * 

The authors wish to express their appreciation for the 
guidance and assistance given them by Professor M. H. Vavra 
of the Department of Aeronautics of the Naval Postgraduate 
School. 
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2. List of Symbols 



A area, in^ 

B number of blades 

b blade height, in 

c p specific heat, BTU/lb^ °R 

c blade chord, in 

D diffusion factor 

DP) driving power 

f r Reynolds number factor 

f^. area multiplier 

2 

g gravitational constant, 32.2 ft/sec 

H total enthalpy, BTU/lb^ 

i incidence angle, degrees 

i Q incidence angle of uncambered blade section, 

degrees 

J conversion factor, 778.2 ft lb/BTU 
Kj^ weight flow blockage factor 

Kj_ correction factor in incidence— angle relation 

K$ correction factor in deviation— angle relation 

1 characteristic length, in 

M Mach number 

m factor in deviation— angle relation at o' = 1 
N revolutions per minute, RPM 

N Re Reynolds number 

n slope factor in incidence— angle relation 
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> /••£ 



n p poly tropic exponent 

P total pressure 

p static pressure 

q dynamic pressure 

q av theoretical average dynamic pressure 

R gas constant/ 53.34 ft Ib/lbjn °R 
r radius/ in 

r.^,^ radius, inner and outer, in 
S entropy, BTU/lb m °R 

s blade spacing, in 

o 

T total temper ature, R 

t> static temperature, °R 

t/c blade maximum— thickness —to -chord ratio 
U peripheral speed, ft/sec 

velocity, ft/sec 
flow rate, lbm/sec 

2 

w equivalent flow rate, in 

X pressure ratio parameter 

x distance from annulus outer wall, in 

z coordinate along axis, in 

« discharge coefficient 

ft air angle, angle between air velocity and axial 

direction, degrees 

X ratio of specific heats 

<£ deviation angle, degrees 

£ density coefficient 

inlet total-pressure -loss coefficient 
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polytropic loss coefficient 
adiabatic efficiency 

pitch angle in model discharge annulus, degrees 
viscosity, lb^/sec ft 
density, slugs/ft 

solidity, ratio of chord to spacing 
dimensionless flow function 
blade camber angle, degrees 

yaw angle in model discharge annulus, degrees 
angular velocity, radi s an,s/sec 
fan total-pressure-loss coefficient 

Subscripts: 

0 station ahead of inlet 

1 station at rotor inlet 

2 station at stator inlet 

2— D low speed two-dimensional cascade 

3 station at stator exit 

10 blade maximum— thickness -to-chord ratio of 10 per- 
cent 

ad adiabatic 

av average 

ax axial 

c compressor 

D design 

d model discharge 



© 

<r 

§ 

<p 

r 

CO 

CO 
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H horizontal 

h hub 

I model entrance 

j summation index 

M model 

m mean 

ma mass averaged 

mer meridional 

min minimum 

n nozzle 

P plenum 

r radial direction 

ST stage 

t tip 

u upstream of nozzle 

V vertical 

z axial direction 

© tangential direction 

Superscripts : 

' relative to rotor 

b exponent in deviation— angle relation 
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3. Design Specifications 

The fan is to deliver air at X = 1.3947 and 
R = 53.34 ft lb/lbro °R. Design point conditions are as 
follows: 

a) Speed = 85,000 RPM 

b) Plow rate = 83.5 lbm/min 

c) Driving power = 23.916 HP 

d) Conditions ahead of fan inlet-scroll-housing 
are: 

1) uniform flow 

2) T q = 7 14°R 

3) P Q = 13.656 psia 

e) Required total pressure ratio -between stator 

discharge and the entrance plane of the 
scroll— housing is P 3 /P Q = l»b78 

The inlet— scroll— housing has a rectangular entrance 
(5.65 in. by 7.82 in.) and a flange 3.125 in. from the cen- 
ter line of the fan. The distance of the center plane of 
the rotor from the back plate cannot exceed 2.1 in. because 
of the critical speed of the shaft. The hub diameter on 
the upstream side of the rotor cannot be less than 1.82 in. 
to accommodate the bearings. The discharge pipe of the 
fan has an ID of 3.75 inches. For manufacturing reasons, 
and to facilitate assembly, the outer diameters of the blad- 
ing should be constant. The fan is to be single-stage and 
without inlet guide vanes. 
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4. Equipment and Instrumentation 

Air supply for the test facility is from a Carrier cen- 
trifugal compressor driven by a 300 HP electric motor. The 
compressor has a rated output of 4 lb/sec at a pressure 
ratio of 2. 

A general schematic of the test installation is shown 
in Figure 1. Hie major portion of the test equipment is 
shown in Figure 2. 

The duct from the plenum to the model is 4 ft. long and 
has 11.3 in. by 15.64 in. inside dimensions. It has two 
3 in. long pieces of honeycomb material to provide uniform 
flow at the entrance to the model. 

Flow rate through the model is controlled with the by- 
pass valves. The primary by— pass provides coarse adjustment 
of flow rate. The secondary by— pass, located in the test 
room, is used for fine flow rate adjustments. 

Flow rate is measured by a VDI standard flow nozzle, in 
the 4 in. inlet pipe, upstream of the plenum. The static 
pressure drop across the nozzle is measured on a 20 in. 
mercury micromanometer. Static pressure upstream of the 
nozzle is measured on a 92 in. water manometer. 

The plenum temperature and the temperature upstream of 
the nozzle are measured with iron-cons tantan thermocouples^ 

Total pressure at the model entrance is measured by two 
Kiel probes mounted horizontally and vertically, and read 
on a 250 cm. water manometer. 
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Measurements in the model discharge annulus are made 
with a United Sensor three-dimensional five hole pressure 
probe. Pressure measurements are read on a 120 cm. water 
manometer. Yaw angle, Y , is read from an angle scale on 
the probe support. 
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5. Inlet Analysis 

Velocity distribution ahead of the rotor and pressure 
loss through the inlet must be known in order to start the 
fan design. Shis determination was made by testing a double- 
scale model of a proposed inlet design. The model casing is 
shown in Figures 3 and 4. The bellmouth for the model is 
shown in Figure 5. i‘The bellmouth hub is shown in Figures 6 
and 7. 

In order to generalize the test results, all flow rates 
were reduced to an equivalent flow rate 




where p and T are the conditions at any point. 

Similarity Conditions 

The ability to apply data from the testing of one size 
machine to another is based on the concept of geometric and 
flow similarity. Geometric similarity was achieved by 
changing all dimensions by the same ratio. Flow similarity 
is defined by the similarity parameters of Mach number and 
Reynolds number. The matching of these two parameters 
assures similar flow. 

Using the definitions of Mach number and the speed of 
sound, the Mach number, in terms of equivalent flow rate, 
becomes 

* 

w 

M = — 7 =- 

Afr 
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Noting that the area of the model is four times the area of 
the design and that is the same, the equivalent flow 

rate of the model, for equal Mach numbers, is 

w (Model) = 4w (Design) 



The Reynolds number is defined by 

g V 1 



%e = 



M 



and in terms of w it becomes 



= 



4 * 

w p 



Re /* VgRT A 



Noting that the area of the model is four times the area of 
the design and that the characteristic length, 1, of the 
model is twice that of the design, the equivalent flow rate 
of the model, for equal Reynolds numbers, is 



w (Model) 



2 P® . v * (Design) 

A<D P m V 



Using the specification conditions for the design, and aver- 
age conditions for the model of 

p M = 14 * 7 psia 

t M = 550° R 

the equivalent flow rate of the model, for equal Reynolds 

numbers, becomes 

w* (Model) = 1.344 w* (Design) 

The equivalent flow rate for the design, calculated from 

2 

the specifications, is 3.505 in. . For Mach number simi- 
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larity this requires a flow rate of 6.8 lb/sec in the model. 
For Reynolds number similarity a model flow rate of 2.29 
lb/sec is required. 

Since the flow through the inlet is in the low subsonic 
range, the effects of Mach number are small and were neg- 
lected* Hie model testing was done at a flow rate which 
gave the Reynolds number similarity condition. 

Test Procedure 

" The flow rate was set and held constant by maintaining a 
constant static pressure upstream of the flowi measuring noz- 
zle. Readings were taken at intervals during each traverse 

to determine the average flow rate and total pressure at the 
model entrance. 

The five hole pressure probe was traversed across the 
annulus height at the model discharge. Measurements were 

taken at each 10% point of the annulus height from the outer 
radius. Additional measurements were taken at- 3, 6, 94, and 
97 percent of the annulus height to get an indication df the 
boundary layer formation. Data was obtained to determine 
total and static pressure, and yaw and pitch angle. Tra- 
verses were taken at eight stations, spaced every 45 degrees 
around the annulus, as shown in Figure 8. 

Data Reduction 

The reduced data for a typical model test is given in 
Appendix A. Due to the volume of data obtained only that 
for Modification 7 (Mod. 7) is given. 
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Equivalent flow rate at the nozzle was determined by the 
method presented by Kelly [ 1 ] . The relevant equation 
from [l ] is 

w n* - * £ f t £ r *n ^ (5-D 

Equivalent flow rate at the model discharge is related 
to the equivalent flow rate at the nozzle by 




where p av is the average static pressure at the discharge. 
Since the discharge velocities were small, the discharge 
static temperature was assumed equal to the total tempera- 
ture, which in turn was considered to be the measured ple- 
num temperature. 

The average static pressures at each of the eight sta- 
tions were determined from the pressure probe data. The 
data from the probe traverses was used in conjtinction with 
the probe calibration chart (Fig. 9) to determine the 
following local flow conditions: the total and static 

pressures, yaw and pitch angles, and the local total- 
pressure— loss coefficient. The parameter used to enter 



Figure 9 is ; > ~ g» 



The further procedures are indicated 



in Appendix A. 

The yaw angle r was measured directly from a vernier 
scale protractor mounted on the probe support. Positive yaw 
was defined as that producing clockwise whirl if looking 
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upstream into the discharge annulus. The pitch angle 
was determined from Figure 9 and defined as positive for 
radial velocity component pointing outward from the dis- 
charge centerline. 

In order to give an indication of loss characteristics 
for various points, stations, and model modifications, the 
local total-pressure-loss coefficient was defined as 

p - P 

< 5 - 3 > 

' P . » p 



where Pj is the average total pressure at the entrance 
to the model. 

In order to compare flow distributions at the various 
stations, and for different modifications, the dynamic 
pressure ratio q/q av was determined at each data point, 
where q av is defined as the theoretical dynamic pressure 
that should occur at the discharge annulus for uniform flow 
and a given discharge equivalent flow rate. That is, 

q = P - p (5.4) 



and 



i-av 






£av 

2 



(5.5) 



Derivation of equation (5.5) was based on the following 
relations: 

continuity equation: w = 0 av gA d V av (5.6) 



equation of state: ^ av = 



'av 



g R t 



(5.7) 



av 
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discharge equivalent 
flow rate: 



w \Tt 



w. 



av 



'av 



4T 



(5.8) 



definition of dynamic pressure: 
(for incompressible flow) 



nq 



av 



>_av v - 
2 av 



(5.9) 



Substitution of V av from equation (5.6) into 
equation (5.9) gives 



„ - 2 
av w 



.2 

w 



aV ^ ?av 2 9 2 *d 3 ' 2 ? av 9 2 A d 2 



or, from equation (5.7) 



w 



Rt 



av 



av 0 2 

2 g A d 



P 



av 



Substituting w from equation (5.8) gives 



L av 



= (w H ) 



*,2 *av 



'av 



g 

R 



R t 



av 

T 



2 g A d Pav 



therefore 



‘av 



■ (?) 



av 

2 



The dynamic pressure ratio then becomes 



JL = 



l av 



w. 



* 



av 

2 



(5. 10) 



The local velocity parameter Vax/ V av was cons idered use- 
ful for noting the effectiveness of the various modifications 
made to the model in an attempt to attain axisymmetiric j flow. 
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For a local velocity V with components in the radial, 
axial, and tangential directions the axial velocity is 
V ax = V cos © cos V • Therefore, for constant density 



Several modifications were made to the model in an 
attempt to attain the same flow conditions at each of the 
eight stations shown in Figure 8. Initial testing indicated 
that deficient flow existed at stations 1, 2, and 8. Sta- 
tions 4, 5, and 6 had excessive flow, while stations 3 and 7 

were receiving the desired amount of flow. Yaw angles were 

o o 

excessive at stations 2, 3, 7, and 8 (about 15 -20 ). Pitch 
angles were highest at stations 4, 5, and 6 (about 8°-10°) . 

In an attempt to correct the flow pattern and achieve 
axisymmetric flow in the model discharge, seven internal 
modifications were made to the model. These modifications 
were tested individually except as noted below. The follow- 
ing is a description of each modification with a discussion 
of the reason for the change and a summary of the results 
obtained. For clarity, the pictures of each modification 
(except Mod. 7) are shown in the relevant figures without 
the hub and bellmouth assembled. 




cos © cos T 



(5.11) 



Modifications 
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Mod. 1 - A flow splitter with a 2 in. radius of curva- 



ture was placed at the top of the model casing (station 1) 
as shown in Figure 10. This change was made for two reasons. 
The original splitter was unavoidably damaged during the 
process of boring a probe hole for stations 1 and 5. It was 
also considered appropriate for the splitter to have some 
curvature so as to provide a smoother change of flow direc- 
tion at the top of the casing. The flow distribution 
(Fig. 11) remained unchanged and the local loss distribution 
(Fig. 12) was almost identical to the preliminary test re- 
sults. For these reasons and due to the fact that this 
modification was used in all subsequent testing/ the Mod. 1 
results were used as a basis for comparison with other 
modifications . 

Mod. 2 — A flow splitter was placed at the bottom of 
the casing (station 5) as shown in Figure 13. It extended 
from the front wall of the casing to the beginning of the 
curvature of the bellmouth. This change was made to provide 
a smoother change of flow path and to direct more flow to 
the top of the casing. The only effect seen in this test 
was an increase in the pressure loss. The flow distribution 
was not altered. Therefor©/ these results are not shown. 

Mod. 3 — A full width flat guide vane was placed in the 
model entrance at the end of the flow splitter of Mod. 2. 

The guide vane alone is shown in Figure 14. It was set at 

zero angle of attack to the flow. Hie trailing edge has a 
circular cutout so that it wraps around a portion of the 
bellmouth. The maximum chord is 4.25 in. at the wall. 
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This guide vane was designed and placed so as to block 
a portion of the flow going to stations 4 , 5, and 6 and to 
force it to the upper stations 1, 2, and 8, where the flow 
was deficient. It was considered to be a gross over- 
correction as it was placed very close to the bellmouth 
lip and extended completely across the model entrance width. 
That is, an over-correction was made in hopes of later re- 
ducing the correction to some optimum. 

The results for this modification showed that high pitch 
angles and high losses existed at station 5. The flow was 
deficient at this station near the outer radius but actual- 
ly increased in magnitude near the hub (inner radius). The 
flow conditions at station 1 were unchanged from those of 
Mod. 1. Therefore, these results are not shown. 

Mod. 4 — Two guide vanes were placed at the model en- 
trance as shown in Figure 15. The vanes are flat and have 
a 1.5 in. chord. They were set at zero angle to the flow , 
and 2 in. from the plane of symmetry of the casing. This 
modification was designed to split the flow into three sec- 
tions. The middle section (between the vanes) would go to 
the lower stations 4 , 5, and 6, and the two outer sections 
would go to the side "and upper stations. The chord length 
was selected so as not to interfere with the flow close to 
the bellmouth. 

The supports for these vanes were designed such that each 
vane could be adjusted over a 2 in. slotted-track and set at 
any angle of attack. For this modification a neutral posi- 
tion was used as a preliminary choice of vane positioning. 
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The results for this test showed no change in flow distribu- 
tion or losses from that of Mod. 1 and are therefore not 
presented. 

Mod. 5 - The guide vanes used in Mod. 4 were reset at a 
25 degree angle to the flow with the leading edge 0.6 in. 
from the plane of symmetry as shown in Figure 16. The re- 
positioning of the guide vanes in this manner was considered 
to be a large correction in a further attempt to reduce the 
flow going to the lower stations and to direct it to the 
upper stations. 

For the purpose of reducing experimental testing time 
such that other modifications could be tested, preliminary 
data for this modification was taken only for five points 
(10, 30, 50, 70, and 90 percent of the annulus height) at 
each of stations 1, 3, 5, and 7. If the results for this 
data indicated significant improvement then a complete 
annulus traverse would have been carried out. 

Comparison to Mod. 1 can be made at the same percentage 
values of the annulus height. Results for Mod. 5 are shown 
in Figures 17 and 18 for axial velocity and loss distribu- 
tions, respectively. The flow at station 1 was increased 
slightly. The flow at station 5 was more uniform but of 
about the same amount as before. Changes at station 3 and 
7 were small. Figure 18 shows that the local losses in- 
creased at station 5 and changed only slightly at stations 1, 
3 , and 7 . 
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Mod. 6 - Two guide vanes were placed at the model en- 
trance as shown in Figure 19. The vanes are "S" shaped with 
a 5.35 in. chord. The leading edge was set at zero angle 
to the flow and was Q.75 in. from the plane of symmetry. 

From the leading edge they have a radius of curvature of 2 in. 
for 2.25 in., then a straight portion for 0.8 in., then an 
opposite radius of curvature of 3.5 in. for 2.55 in. to the 
trailing edge. The vanes are tangent to the bellmouth at 
the trailing edge. 

This modification was a gross correction to divide the 
flow into three sections and force sufficient flow to the 
upper stations. Although considered impractical this change 
was made and tested in the hope that an optimum and more 
practical configuration could be established. 

The Mod. 6 Velocity Distribution (Fig. 20) shows that in 
comparison with Mod. 1 the flow was increased at all sta- 
tions except 4 and 6 which were very deficient. Stations 1, 

2, and 8 now had excessive flow, and the flow at station 5 
became even more excessive. Several stations had less uni- 
form axial velocity distributions than in Mod. 1. Stations 4 
and 6 experienced very high losses as shown in the Mod. 6 
Loss Distribution (Figi*.r:>21.}^‘ t 

Mod. 7 - A groove 1 in. deep and 2.28 in. long was cut 
in the bellmouth outer wall as shown in Figures 5 and 22. 

The purpose of this change was to provide more area for flow 

to the upper stations and to improve the bellmouth shape. 

The assembled configuration is shown in Figure 23. 
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Figures 24 and 25 show the axial velocity and loss dis- 
tributions, respectively. The results in general are very 
similar to those for Mod. 1. However, the station that 
changed the most was station 1, and this was as an improved 
(increased) flow distribution. Stations 2 and 8 showed 
increased losses near the outer radius. Station 1 experi- 
enced somewhat lower and more uniform losses. The other 
stations indicated little or no change over Mod. 1. The 
pitch and yaw angles were reduced slightly, but fairly sig- 
nificant yaw angles (15°) still existed at stations 2, 3, 7, 
and 8. 



Predicted Conditions 
Ahead of Fan Rotor 

The final part of the inlet analysis was the prediction 
of the flow conditions ahead of the rotor for the fan design. 
This prediction was based on the non-dimensional velocity 
ratio and the total-pressure-loss coefficient from the model 
tests, and on the fan-inlet entrance conditions given in the 
specifications. 

Final Model Analysis. - The model test data indicated 
that the Mod. 1 and Mod. 7 configurations were similar. For 
better comparison the data for these two conf igurations was 
circumferentially averaged for © , Y • 7 ' and V ax /V 

av* 

It was decided that Mod. . 7 gave slightly better flow distri- 
bution between the outer and inner walls, and that the 
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averaged yaw and pitch angles were somewhat smaller. The 
loss characteristics of the two configurations were nearly 
identical. Therefore, the Mod. 7 data was chosen for use 
in the prediction of rotor entrance conditions. 

The circumferentially averaged velocity distribution, as 
shown in Figure 26, was somewhat in error as the area under 
this curve was about three percent greater than unity 
(1.0316), although theoretically this area must equal unity; 
that is, 



Conditions that were considered likely causes of this dis- 
crepancy are: decreased discharge area due to probe block- 

age, general instrumentation errors, improper estimation of 
the velocity profile near the walls, increased probe in- 
accuracies in the boundary layer, and the assumptions made 
in the calculations (such as t^ = T^ , ^ = constant, etc.). 

To correct this apparent error the Mod. 7 velocity dis- 
tribution was normalized such that equation (5.12) was satis- 
fied. This was done by dividing each ordinate value of 
V ax /V av by the area under the curve. Both distributions, 
Mod. 7 and corrected Mod. 7, are shown in Figure 26. 

The mass-averaged total-pressure-loss coefficient for the 
model was determined as follows 



1.0 




(5.12) 




(5.13) 
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is the average total pressure at the model en- 



where P^ 

trance, and P^ and p nia are the mass-averaged total and 
static discharge pressures, respectively. Thus, 



The static pressure p at discharge was mass-averaged 
in the same manner. The integration was carried out radi- 
ally and then circumferentially. Using equation (5.13) it 



Design Prediction . - The conditions ahead of the rotor 
were determined by the method of Vavra [ 2 ] for an adia- 
batic expansion process with friction and varying cross- 
section. 

The fan-inlet loss coefficient was assumed to be equal 
to the mass-averaged value for the model. 

The dimensionless flow function was determined from 

the relation in £ 2 ] 



P, 



> — S dw P 

ma Jdw 



which may be written as 



P 





V av dr d 0 



For assumed constant density, 





was determined that ^ ma = 0-0484. 
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w 



1 



(5.14) 



$ 



P, 



A n 




R 



where w, T 0 , and P Q are given in the design specifica- 
tions. The area ahead of the rotor is dependent upon 

the design choice since dimensional restrictions are given 
in the specifications. The double-sized model was based on 
the minimum hub distance and a casing diameter of 3.5 inches. 
Thus, the value of Aj_ was prescribed. 

With and as known parameters, Table Cl of £ 2 ] 

was used to find the pressure ratio P Q /p^ . This procedure 
was based on the simplifying assumptions that 



So - 1.4 

and 

* > D = 0.05 

The following equations from [2] were then used to 
determine the velocity : 




where rip is defined as the poly tropic exponent, ?e as the 
polytropic loss coefficient, and X as a pressure relation 
used as a simplification of the equations. 
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The velocity represents the required average velo- 

city ahead of the rotor for a non-uniform velocity distri- 



bution, or the actual velocity for uniform flow that is 
entirely axial. 

Since = T Q for an adiabatic process, the static 

temperature was determined by 



t 



1 




2g J Cp 



The total pressure P-^ was determined by 



where 



P, 





P 



o 





1 



The design weight flow was used as a check on these pre- 
dicted conditions by 




A 1 V 1 



The absolute Mach number is 

% 

M 1 = V ]/ 9. R fc l) 

Given below is a summary of the properties at the rotor 
inlet: 

P x = 13.534 psia T x = 714° R V x = 627.6 ft/sec 

p-^ = 11.483 psia t-^ = 681.55° R = 0.490 
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6. Aerodynamic Fan Design 



The basic equations that are applicable to the design of 
an axial-flow compressor or fan are developed in £ 3 ^ . 

These are formulated from the laws of conservation of matter, 
momentum, and energy, along with the thermodynamic equation 
of state of the working fluid. 

The simplifying assumptions that are made in 3 } are 
as follows: 

1) The general flow equations are applied only to 
compute flow distributions between blade-rows 
where blade forces are nonexistent. 

2) The flow is assumed to be steady and axially 
symmetric. 

3) Viscosity effects between blade rows are neg- 
lected in the flow equations. However, these 
effects are considered by use of experimentally 
determined blade-element performance data. 
Empirical corrections are also made for wall 
boundary-layer effects on required flow area. 

4) Heat transfer is neglected. 

The two principal phases of the aerodynamic fan design 

are: 

1) Determination of blade-row velocity-diagrams 
for design point operation. 

2) Selection of blading to satisfy these velocity- 
diagrams. 
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Determination of Velocity -Diagrams 



General . - The determination of the various air veloci- 
ties and flow angles from hub to tip at the inlet and outlet 
of the rotor and stator was made based on the radial equili- 
brium equation [ 3 J , which is 



JgC p ^ _ Jgt If + 






V. 



r 



^v 2 

v z 3T 



- V, 



^v r 

J7 



( 6 . 1 ) 



This equation results from the preceding simplifying assump 
tions, and was the theoretical starting point for the fan 
design. 

It was assumed that the meridional velocity 

I 2 2 ~ 

V mer = v v z + v r is 9i ven by the axial velocity V z , 
and that ^V r /<Jz = 0. Such a condition has been success- 
fully used in low-aspect-ratio and lightly loaded blade- 
row designs. Therefore, this assumption was considered 
appropriate for the fan design. 

It was decided that the design should include constant 
energy addition through the rotor blade-row. That is, the 
change of total enthalpy through the blade-row is constant 
along a streamline. If the total enthalpy is assumed con- 
stant ahead of the rotor (H^ = constant) , then a constant 
total enthalpy will exist behind the rotor. This may be 
expressed as 




The Euler turbine equation may be written as 
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H 2 “ H 1 = Jg t ^ rV ® * 2 “ ^ rV ® ( 6 * 2 ) 

Since V e =0 ahead of the rotor, and Hj - is constant, 

then (rV 0 ) 2 = constant. Thus, the above conditions re- 
quire that free-vortex flow exists behind the rotor blade- 
row. 



Therefore, neglecting radial velocity components, equa- 
tion (6.1) becomes 



V fl 



*V, 

w 



& 



V, 



+ v + t = 0 

2 3ft or 



(6.3) 



The first two terms in this equation may be written as 




or 

V 0 (r V e ) 
~r~ ^ r 



Since rV Q = constant, the sum of the first two terms in 
equation (6.3) equals zero. 

The radial equilibrium equation then reduces to 



a v z Js 

v *ir +t 5r 0 



(6.4) 



The simplest solution of equation (6.4) is the "simple- 
radial-equilibrium solution, " based on the further assump- 
tion that the derivative of S with respect to r is zero. 
This condition has been widely used in compressor design, 
and was considered valid for the single-stage fan for slight 
change in annular area. This assumption is made only at a 
fixed value of z . 
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Therefore, equation (6.4) becomes 



with the solution 



V, 



J v z 

^r 



0 



V z = constant 



(6.5) 



Using this result as a basis, the velocity diagrams and 
flow properties were calculated for two assumed rotor inlet 
conditions; namely, 

a) Design A for the circumferentially-averaged 
axial velocity for the most satisfactory flow 
determined from the inlet model tests. The 
mass-averaged total-pressure-loss coefficient 
for the model was assumed. Total temperature 
and total pressure were assumed to be constant 
ahead of the rotor. Pitch and yaw components 
of velocity were assumed to be zero. 

b) Design B for uniform flow with axial velocity 
components only. The velocity and properties 
ahead of the rotor were based on the same loss 
coefficient as Design A. 

General calculation procedure . - The calculation method 
is one of iteration and compromise. Iteration is necessary 
in order that the continuity equation may be satisfied by 
proper annular dimensions. Compromise is necessary in order 
that blade-loading and Mach number limits are not exceeded. 

Experimental tests of axial-flow compressors show that a 
Mach number of 1.2 is the approximate limit for a transonic 
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design method as used in this paper. Further, the following 
diffusion-factor limits are appropriate for inlet stages: 
for the rotor tip, D less than 0.4; for the rotor hub, b less 
than 0.6; for the stator, D less than 0.6 [ 3 ] . Although 
these criteria are somewhat vague, they remain the best 
indicators of design limitations that are currently avail- 
able. 

The detailed development of the calculation method is 
presented in [ 3 ] «• The iterative procedure is started by 

selecting the radial stations according to percentages of the 
passage depth. The computation scheme and results for velo- 
city calculations are shown ih Tables I and II, for Design A 
and Design B, respectively. Only the final results and the 
final calculation lineups of the iteration procedure are 
presented. 

The velocity leaving the stator was assumed to be entire- 
ly axial, thus providing maximum flow rate for a given stage 
loading. 

Design A. - Table I shows the steps taken to determine 
each parameter. Figure 27 shows a general stage velocity- 
diagram which is indicative of the fan design. 

From the results of the inlet-model analysis, the assumed 
inlet conditions to the rotor were as follows: 

T x = 714° R 

P-j_ = 13.534 psia 

V z l 

> — as shown in Figure 26 

V av 
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V av = 627.6 ft/sec 

t = 1.74 in. (for an assumed tip clearance of 0.010 

in. ) 

r. . = 0.91 in. 
l,h 

The conditions at station 1 are listed as steps 1 through 
8 in Table I. The iteration was started by assuming that 
the hub radius was constant through the rotor (step 9) . Most 
of the steps are clearly indicated under the procedure col- 
umn in Table I, in which Figure 27 was used. However, sev- 
eral steps had to be related to important design assumptions 
and will be discussed in detail. 

V (step 11) was determined by use of the Euler rela- 
i y / 2. 

tion, equation (6.2), and the definition of DP, the driving 
power. That is, 

DP = w(H 2 - 

or DP = ^ (rV e ) 2 

Jg 2 

Since oj = 7TN/30, the tangential velocity behind the 
rotor may be found as a function of r 2 . 

It has been pointed out in £ 3 ] that large axial velocity 
reductions across any blade-element are not conducive to high 
efficiency. In addition, the available experimental data for 
transonic rotors shows that the range of V_ 0 to V . ratios 
across the rotor are as follows: for the tip, 0.85 to 1.05; 
for the mean radius, 0.90 to 1.15; for the hub, 0.95 to 1.50. 
Therefore, in step 14 the axial velocity through the rotor 
at the mean radius was assumed to be constant. This assump- 
tion causes the axial velocity changes near the tip and hub 
regions to be typical of the experimental data. 
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Steps 24 through 26 involved a compromise in which sever- 
al important design conditions were considered. These were 
as follows: the critical blade-chord Reynolds number, fric- 

tional losses, shock losses, blade-loading, manufacturing 
restrictions, and the practicality of rotor diameter measure- 
ment. 

The critical blade-chord Reynolds number may be defined 
as the limiting value of Reynolds number below which a large 
rise in loss occurs. Above this value the losses vary only 
slightly or not at all. The available experimental data are 
generally for flow conditions at or above; this critical 
Reynolds number. Therefore, rotor blade design that is 
based largely on test data should also be for similar flow 
conditions. Based on the relative inlet velocity to a 
blade-element, Vavra [4 } indicates that this limiting 
Reynolds number is about 2 x 10^. 

Too high a solidity is associated with increased fric- 
tional losses. Too low a solidity results in higher blade 
loading, and increased shock losses near the tip. Recom- 
mended rotor tip solidity for transonic rotors in regard to 
shock losses is 1.0 [ 3 ] . 

For manufacturing reasons the minimum blade spacing at 
the hub is about 0.25 inches. 

An even number of rotor blades allows for easier measure- 
ment of the tip diameter. 

Therefore, in consideration of the above facts, the mini- 

5 

mum blade-chord Reynolds number was set at 2 x 10 for the 
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rotor. 

mined. 



The corresponding minimum blade-chord was deter- 
for the known conditions ahead of the rotor, by 



%e “ 



v i c 



fi 






( 6 . 6 ) 



where was found in [ 5 ] for known values of t-^ . 

Several combinations of chord distributions and blade 
numbers were investigated in order to determine a satisfac- 
tory solidity range. Pertinent equations used were 



2 TTr 

s = — 



and 



0 * 



c 

s 



It was considered important to have a tip solidity fairly 
close to 1.0 as recommended in [ 3 3 / but also essential to 
establish a solidity range similar to that for the experi- 
mental data for transonic rotors. 

As a result, 12 rotor blades with uniformly varying 
chord were selected such that the Reynolds number was at or 



above 2 x 


10 5 . A 


summary 


of this analysis. 


which was 


based 


on no area 


change 


through 


the rotor. 


is as 


follows: 




Radial position 


V 1 


c min ° 


! t J 

, 

'9 ; 


s 


<r 


TIP 




1450.1 


0.536 


0.750 


0.911 


0.823 


. 10 % 




1395.5 


0.557 


0.760 


0.868 


0.876 


50 % 




1169. 1 


0.665 


0.800 


0.694 


1. 153 


90% 




960.4 


0.809 


0.840 


0.520 


1. 615 


HOB 




914.0 


0.850 


0.850 


0.476 


1.786 



From this summary the linear chord distribution was ex- 
pressed as a function of radius for those iterations to fol- 
low. Thus , 
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c 



(6.7) 



= c 



l,h 



c 

l,h 



r 



l,t 




(r-r 



l,h 



) 



It should be noted that the above chord-solidity summary 
refers only to the first iteration. With reduced blade 
height at the rotor discharge the values of r-^ and ^ 
will be slightly different for a given percentage of passage 
depth. A streamline will follow a curved path through the 
rotor, depending up^>n annulus constriction. Therefore, the 
solidity must be averaged for the values at rotor entrance 
and exit. This was done in step 27, where o' 1 is referred 
to the chord and spacing at station 1 for a given percentage 
passage depth from the casing, and O ' 2 to the chord and 

spacing at station 2 for the same percentage passage depth. 

/ 

The total-pressure -loss coefficient co was determined 
from the experimental data for a geometrically similar 

double-circular-arc transonic rotor £ 6 3 • This rotor was 

/ 

tested under similar flow conditions. The values of <JJ 
were for minimum-loss incidence angle. This method of 
determining the loss coefficient was considered appropriate 
for the simplified theory used in this paper. 

The relations used in steps 30 to 32 are developed in - 

1 3 ] . 

It was considered probable that the boundary layer in the 
fan rotor would have significant blockage effects, especially 
due to the small blade heights. Therefore, in step 34 a 
blockage factor of 0.96 was assumed. The continuity equation 
was checked by graphical integration. 
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The iteration was repeated for a new hub radius at sta-» 
tion 2, re-starting at step 9 until continuity was satisfied. 

The final rotor outlet conditions were used as stator 
inlet conditions. A similar process was followed for the 
stator design. However, the stator blade solidity was deter- 
mined in a somewhat different manner than for the rotor, as 
excessive solidities resulted for the criterion that minimum 
Reynolds number be equal to 2 x 10^. 

A solidity distribution was sought which would be simi- 
lar to the available experimental data for the NACA-65 (A-^,) - 
series airfoil. (The specific reasons for choosing this 
airfoil shape will be given later) . The extensive test data 
for this airfoil includes solidities up to about 1.6, but 
extrapolated data for blade solidities up to 2.0 are given 
in £ 3 ] . 

Several combinations of blade-chord and number of blades 
were investigated. As a result, the solidity at the 90% 
position (near the hub) was set at 1.55, and the solidity 
at the 10% position (near the tip) was set at 1.0. Blade 
chord was chosen as constant at 0.805 inches, and the number 
of blades was set at 13. 

Step 42 indicates the resulting Reynolds numbers. These, 
values were deemed acceptable for two reasons Q 3 3 • First, 

c 

the relative difference between Reynolds numbers of 1.5 x 10 
ana 2.0 x 10 5 is small for the NACA-65 (A 10 ) - series airfoil 
under moderate blade loadings. Secondly, lightly loaded 
blades (such as the fan design of this paper) have only a 
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gradual increase in loss below the critical Reynolds number, 
rather than the more typical steep rise in loss for normal 
blade loadings. The diffusion factor values listed in step 
43 indicate that the fan stator is rather lightly loaded as 
compared to the suggested design limit for stators of D less 
than 0.6 . 

The total-pressure-loss parameter (step 44) was deter- 
mined from £3] for known diffusion factors. 

The continuity equation was checked using an assumed 
blockage factor of 0.94. Iteration was found to be unneces- 
sary. 

The stage total pressure ratios and stage adiabatic effi- 
ciencies for the three radial positions were determined. 

These values were mass averaged as shown in Table I. 

Finally, the total pressure ratio over the inlet-fan 
combination was determined for comparison with the design 
specification. 

Design B. - Table II shows the procedural steps taken 
to determine each parameter. 

From the results of the inlet analysis for uniform flow, 
the conditions ahead of the rotor were taken as: 

V 1 = v z,l = V av of DESIGN A = 627.6 ft/sec (uniform) 

t 1 = 681.5° R 

Total properties and rotor radii the same as Design A. 

The detailed steps and design assumptions were similar to 
Design A. Iteration and mass-averaging were likewise similar- 
ly done. 
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Selection of Blading 



General , - The blade-sections at each of the three 
radial positions were selected in order to satisfy the de- 
sign velocity-diagrams. 

Double-circular-arc blading was chosen for the rotor be- 
cause this relatively simple airfoil shape has been used 
effectively in transonic compressors. 

The NACA-65 (A^q) -series airfoil was considered to be 
appropriate for the stator blading. The reasons for this 
choice were the magnitudes of the design inlet Mach numbers 
for the stator and the availability of extensive experiment- 
al data for this airfoil. 

The parameters that describe blade-element flow are 
total-pressure loss, incidence angle, and deviation angle. 
The total-pressure-loss coefficient for minimum loss inci- 
dence has already been established in the determination of 
the velocity-diagrams. Therefore, the design incidence 
angle is for minimum loss. 

Camber angle, incidence angle, and deviation angle are 
used to define the blade camber, air approach, and air 
leaving directions, respectively. These angles are shown in 
Fig. 28 for a typical blade-element. 

Rotor. - A rotor-blade selection procedure is presented 
in C 3 3 for the double-circular-arc airfoil. This method 
is based on a correlation of experimental test results with 
both annular and rectilinear cascades. Empirical correction 
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